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Abstract 

The standard Frohlich model of the interaction of band states with phonons is 
applied to the case of deuterons solved in metals. A simple variational solution 
of this model predicts at least one highly coherent state of the deuteron-phonon 
system having a total energy significantly lower than that of uncorrelated interstitial 
deuterons and only thermal phonons, and that this phase may actually be exhibited 
at room temperature and high levels of deuteration. 
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1 Introduction 



A recent attempt by Kalman and Keszthelyi [1] to explain the anomalously 
large fusion cross-section observed when low energy deuterons traverse Pd foils 
at low energy was found insufficient to account for the full magnitude of the 
reported effect. The present work was motivated by the observation that their 
model is based upon an a priori transformation of the phonon-deuteron vertex 
of the Frohlich Hamiltonian into an effective deuteron-deuteron interaction, 
and is therefore unable to take account of any collective effects that might 
arise from coherent phonon excitations. 
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2 Model 



The essential physics of hydrogen nucleii solved in metal lattices can be cap- 
tured by a single-band Frohlich model because the lowest states available form 
a band, albeit a very narrow one [2], that is separated from the excited states 
by a constant band gap that is well in excess of both kT and ujd- 
The Hamiltonian for a variable number of deuterons solved on a lattice in 
exchange equilibrium with the gas phase is 

H = e k c k ] c k + M q (b q + bl q )c\ +q c k + U q c\ i+q c\ 2 „ q c k2 c kl + u q b\b q (1) 



e k is the band energy of a deuteron of crystal momentum k and is negligible 
in comparison with the other terms, c k is the raising operator for a deuteron 
of wavevector k, b q is the raising operator for a phonon of wavevector q. We 
have neglected spin and band indices for the sake of clarity. 
According to the Thomas-Fermi model [3] of electronic screening in metals 



q ~ q 2 + 4ne 2 N{e F ) 



The bosonic character of the deuterons is preserved in the metallic environ- 
ment because the accompanying electronic cloud is made of very many small 
contributions with zero average spin. 

Consistency with the classical model of the dielectric response of the lattice, 
and consideration of the fact that the retarded interactions will be subject to 
exactly the same electronic screening response as the instantaneous Coulomb 
interaction implies that 

\M q \ 2 « u q U q (3) 



The phonon vertex interaction is clearly strongest when the deuterons equally 
populate two distinct states separated by the momentum corresponding to 
half a reciprocal lattice vector, G/2 . 

We denote the state having N d d spin- up deuterons with k = 0, an equal 
number of spin-down deuterons with k = G/2 and N p phonons with q = 
G/2 by \N dd ,0,0,0,0, N dd , N p ) . The subspace spanned by this state together 
with the phonon-coupled state \N dd — 1, 0, 0, 0, 1, N dd , N p — 1) has the explicit 
Hamiltonian: 

1 2N 2 dd U + N p u D jN~ p N dd M G/2 
{ \J~N~ P N dd M G /2 (2N 2 d -2N dd + 2)U + (N p -l)cu D 
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where ujd is the Debye energy of the acoustic phonon branch. 



3 Numerical estimate 

With the value of N(e F ) = 0.47 states/(eV Pd) for PdH reported in [3] we 
have U ~ 0.8eV and Uq/2 ~ 0.4eV. For the acoustic branch in Pd, ujd ~ 
0.03eV. Substituting the dimensionless parameters i) mean deuterons-per-site 
A = 2Ndd/N c and ii) mean squared phonons-per-site v 2 = N p /N c - and using 
(2) we obtain in the limit of large N c , in eV units: 



In view of the variational principle, this calculation represents an upper bound 
for the lowest energy e_. If the subspace is expanded to include all even val- 
ues of N p , for example, the splitting caused by the phonon-interaction term 
is found to be about twice as large as that given by (4) above and the corre- 
sponding site energies (in eV) are: 



With a chemically achievable A = 1.0, the lowest value of e_ ^ — 2eV is 
achieved when v m 6. This represents a significant kT) gain in energy 
relative to (normal) uncorrelated interstitial deuteron states without coherent 
phonon excitation. It should also be noted that if the metal is in exchange 
equilibrium with the D 2 gas phase, (6) predicts a spontaneous take-up of ad- 
ditional deuterons once a certain threshold value of A has been reached. 



4 Conclusion 

We conclude that coherent phonon-deuteron interactions may dominate the 
Coulomb repulsion energy in certain deuterated metal regimes, even allowing 
for the phonon energy invested in the highly-excited phonon mode. We also 
believe that the principle outlined here, which is closely related to the well- 
known Peierls instability, may be applicable to other situations involving high 
concentrations of light charged bosonic interstitials in rigid lattices. 





e± = 0.4A 2 + 0.03z/ 2 ± 0.36z/A 



(6) 
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